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Plasmonic quasicrystals (PlQCs), by integrating the properties of quasicrystals (rotational symmetry and 
long range ordering but lack translational symmetry) and surface plasmon polariton mediated effects, offer 
several advantages over plasmonic crystals (PlCs). For example, in PlQCs one could have broadband, 
polarization independent response. However, large area patterning by electron beam lithography requires 
precise lattice coordinates as well as a practical way to design the structures for specific spectral response. We 
demonstrate design and fabrication of large area quasicrystal air hole patterns of n/5 symmetry in metal film 
in which broadband, polarization and launch angle independent transmission enhancement is observed. We 
demonstrate bi-grating quasicrystals to show that designable transmission response is possible over visible 
to near infrared wavelength regions with about 15 times enhancement. These would be useful in many 
applications like energy harvesting, nonlinear optics and quantum plasmonics. 



Sub -wavelength patterns in the form of grooves, hole arrays, or corrugations on metal-dielectric interfaces 
help optically excite surface plasmon polaritons (SPPs) at the interface 1 . SPPs enhance transmission as well 
as the local field and have wide ranging applications in nanophotonics 2 " 12 . Plasmonic crystals are periodic 
metal- dielectric structures realized for a specific application in which the SPPs are optically excited by periodic 
structures 1 . The translation symmetry of these crystals in a given direction results in discrete SPP excitation. For 
example, the longitudinal SPP modes are excited by 1 -dimensional patterns (gratings) by the transverse magnetic 
(TM) field though under specific orientations transverse electric field may excite SPPs 13 . In the case of 2-d patterns 
with circular constituents (holes or discs), while for square lattice the TE and TM polarizations are equivalent due 
to n/2 rotational symmetry, for non-square lattices and lattices with asymmetric shaped holes, unless specifically 
designed, the response will be different for different polarizations for non-normal incidence. In general, the 
projection of momentum onto the surface is given by k x 2 + k y 2 = k SPP 2 , where k x = kosin(0)cos((p) ± i27i/a x and 
k y = kosin(0)sin((p) ± j27r/a y in which i, j are integers, a x , a y are lattice constants in the in-plane directions, 0 is the 
angle of incidence and (p is the azimuthal angle. It would be interesting if the dependence of the discrete SPP mode 
excitation on the orientation of the sample or the azimuthal angle cp and thus polarization of light can be 
overcome 1516 . Unlike crystals, quasicrystals have only long range ordering and rotational symmetry 14 . In quasi 
crystal patterns, discrete SPP modes 17 , extra- ordinary transmission in different wavelength regions 18 " 24 , higher- 
harmonic generation 25,26 , focusing and directivity of emission 27,28 have been demonstrated. In addition to these 
designable k-space and broadband local field and transmission enhancement in plasmonic quasicrystals (PlQCs) 
would be useful. For example, for light harvesting, the top contact of the photovoltaic device can be patterned with 
quasicrystal structure, in nonlinear optics for harmonic generation though the metal surface optical nonlinearity 
is weak (x (2) ) 29 > in quantum plasmonics for high quality factor cavities 30 . Though in random patterns one can 
achieve broadband response 31 , it is in practice not easy to design them for specific k-space response. Aperiodic 
structures and their applications to nanophotonics have been reported earlier 32,33 . 

In order to design quasicrystal structures several methods based on tiling, grid, inflation/deflation, projection, 
decoration, cluster model, etc., are employed 32,34 . In most of the previous reports on PlQCs, the dual grid method 
or the Penrose tiling using fat and thin rhombus shaped tiles for designing and ion beam milling for patterning the 
thin metal film were employed 18 " 28 . Electron beam lithography has been used for making gold nanodots arranged 
in aperiodic patterns 31 and Moire nanolithography has been used to demonstrate quasicrystal patterns of high 
rotational symmetry 35 . In optical domain, the effect of size of the pattern on the transmission enhancement was 
studied as well as the origin of transmission enhancement was proposed to be not related to long range order- 
ing 18,19 . Both these studies are on patterns covering an area of 10 X 10 um 2 milled in silver thin films and thus are 
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comparable to or smaller than the propagation length of SPPs. We 
present a method to realize quasicrystal patterns over large area using 
electron beam lithography. The design and fabrication methods can 
be easily extended to any rotational symmetry and Penrose tiles. We 
present in PlQCs with n/5 symmetry, polarization independent, near 
dispersionless broadband transmission enhancement compared to 
unpatterned metal. We demonstrate ultra broadband transmission 
enhancement in PlQCs. 

Results 

In oblique tiling method, a set of parallel lines tilted by n/n angle are 
superimposed by n similar sets, each rotated by n/n with respect to 
the previous set. The intersection points of the n sets of these parallel 
lines give the coordinates of the quasicrystal with n-fold symmetry. 
For generating the soft mask for electron beam lithography, the 
coordinates of the quasicrystal lattice points are needed. For this 
purpose, we modify the design method by replacing the continuous 
lines with dots and thus the set of parallel lines are now a 2-D array of 
dots as shown in Figure la. By successively superimposing n such 2- 
dimensional arrays of dots that are rotated by n/n with respect to the 
previous array, we generate a hole pattern. The spatial points which 
are common to all the n overlays give the coordinates of the QC of 
n/n rotational symmetry. 

There are two main advantages of this design methodology. The 
intrinsic period of the starting 2-d array defines the center of the 
spectral response and thus can be tuned to any wavelength required. 
Thus, the spectral response of the quasicrystal pattern is designable. 
Secondly, this method simplifies the complexity in identifying the 
coordinates of the quasicrystal lattice points from the large set of 
points for generating the CAD for electron beam lithography. 
However, the repeated points need to be removed to avoid multiple 
exposures. To find repetitions in the large array of points, we used 
Quicksort algorithm which is a sort-in-place algorithm that does not 
need additional memory and sorts an array of size N in 0(N X 
log(N)) time 36 . This helps identify repetitions in approximately 
O(N) time which is much faster than the time taken by brute force 
methods (0(N 2 )). 

The coordinates of the quasicrystal pattern are generated over an 
area of ~ 1 mm 2 . By electron beam lithography and dry etching, we 
patterned thin gold film on optically flat quartz (see Methods). 
Fig. lb shows the SEM image of the fabricated quasicrystal. The 
initial rhombus used in designing the structure had 1 mm long sides 
each having holes at a periodic separation (a in Fig. la) of 600 nm. 
That is, in each set, there are about 2.8 million points in 1 mm 2 area 
and with 5 superimpositions there are a total of about 14 million 
points. 



In the following we use PlQCs in which all the holes excepting 
those that are within 10 nm vicinity of a given hole are retained to 
avoid over exposure during e-beam lithography. This sorting results 
in, the coordinates of both QC of n/5 symmetry as well as those 
corresponding to the 5 base lattices, about 13 million holes in 
1 mm 2 area. The average hole diameter obtained was 92 ± 7 nm 
by careful e-beam dosage. The diffraction pattern (Fig. lc) recorded 
using a 635 nm diode laser beam clearly shows n/5 symmetry. The 
observed multiple concentric circles show that the structure displays 
ordering at different levels. Fig. Id shows the calculated Fourier 
transform of the structure shown in Fig. la. 

We performed angle resolved white light transmission measure- 
ments to measure the normalized transmission (ratio of transmission 
through patterned and unpatterned Gold) for TM, TE and unpolar- 
ized light for different 0. Fig. 2(a) shows the measurement geometry 
and Fig. 2(b) compares the spectra for different polarizations for 6 = 
5° and (p = 0°. Fig. 2(c) compares the transmission enhancement 
spectra for different 0 and (p. Transmission spectra measured for 
different 0 are plotted as contour plots for TE and TM polarizations 
in Fig. 2(d) and Fig. 2(e), respectively. Transmission dips correspond 
to SPP absorption and thus the dark contours show dispersion of 
various SPP modes. Similarly, bright regions correspond to trans- 
mission enhancement regions. We observe upto 9 fold transmission 
enhancement with respect to the unpatterned metal tranmission and 
near dispersionless broadband response seen as the featureless bright 
band between 550 nm and 720 nm wavelengths. 

Discussion 

To calculate the SPP resonances and their dispersion, for a periodic 
pattern, there are theoretical models for a given interface 371 . Analytical 
expressions for SPP dispersion for two layers (single interface) and 
three layers (two interfaces) are well known 35 . In order to calculate the 
grating excited SPP modes one can equate the 2-d momentum con- 
servation equation where appropriate k S pp is used depending on the 
layer structure. For example, for 2-d pattern of air holes 38 , 



y yc 0 sin(0) cos(#?) + —^—j + sin(0) sin(^) ± ~^—J = fapp 

where ai and a 2 are the periods in x and y directions, respectively. 
Compared to a periodic crystal, the azimuthal angle (cp) dependence in 
a quasicrystal is negligible due to rotational symmetry. In addition, 
due to lack of translational symmetry, the Bloch periodic terms {mini 
ai and n27t/a 2 ) do not exist in the quasi crystals and are to be replaced 
with the reciprocal vectors which can be obtained from either the 
calculated k- space picture or from measured diffraction patterns. 





Figure 1 | (a) Schematic showing the steps in designing of quasicrystal with tt/5 rotational symmetry (see text for details), (b) SEM image of the 
fabricated quasicrystal pattern, (c) Measured diffraction pattern of the quasi crystal structure for 635 nm light, (d) Calculated Fourier transform of the 
final pattern shown in (a). 
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Figure 2 | (a) Schematic of the measurement geometry. Normalized transmission spectra are shown for different (b) incident polarizations and (c) 9 and 
(|) for TM polarized light. Measured launch angle dependence of the transmission enhancement is shown as contour plot for (d) TM and (e) TE 
polarizations. Top scale bar shows the transmission enhancement factor. Broadband transmission enhancement (red regions between 500 -700 nm 
wavelength region) that is launch angle independent is seen. Some of the dispersive SPP modes are identified from the calculations. 



Thus, the momentum conservation rule for QCs is given by 17 , k\ \ + 
k SPP = F (i) where k\ \ and k SPP are the in-plane components of the wave 
vector of incident light and SPP, respectively and F (i) are the reciprocal 
vectors. From the calculated Fourier transform of the pattern, one can 
extract the F (i) s and thus calculate the k SPP s. We use this method in the 
following to calculate the SPP resonance positions which can, in prin- 
ciple, be used to calculate the expected spectral response. In practice, 
the dense k- space makes it difficult to calculate the spectral response. 
Another practical advantage in calculating the SPP resonances is, if the 
calculated SPP resonances do not match the requirement, we can scale 
the pattern for specific k- response by using the scaling property of 
Fourier space which, for a constant c and real and phase space func- 
tions given by F and f, is given by, j F(ck)e lkx dk= ~f(^~J- 



We use the Fourier transform (FT) shown in Fig. Id to deduce the 
reciprocal vectors (k l5 k n etc.,) corresponding to the successive cir- 
cles. For example, the values are k x = 1.82 |im _1 (red) and k n = 
2.91 |im _1 (green). Each SPP mode may be denoted in terms of a 
tuple (K, B, m, n). K takes the values k l5 k n , B takes the values Q 
(quartz) or A (air) which defines the interface with gold, at which the 
SPP propagates and (m,n) are integers. The calculated positions of 
various modes are marked in the measured SPP dispersion contour 
plots (Fig. 2(d) and Fig. 2(e)). For the calculations, the values of s m 
are taken from Johnson & Christy parameters 39 for gold and s d = 1, 
2.25 for air and quartz, respectively. The resonance corresponding to 
k n is more prominent as can be seen in all the plots. This is consistent 
with the stronger diffraction maxima in Fig. Id for k n compared to k^ 

Fig. 2(b) shows the broadband transmission enhancement for TE, 
TM and unpolarized incident light for 0 = 5° and (p = 0°. 




400 450 500 550 600 650 700 750 800 850 

Wavelength(nm) 

Figure 3 | (a) Measured diffraction pattern over a square lattice with identical hole density shows much sparser and weaker k-space compared to the 
quasicrystal pattern, (b) Shows the transmission enhancement spectrum for square lattice of air holes in gold. 
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Figure 4 | A comparison of the transmission enhancement for two PlQCs 
of 1 um period with differing inhomogeneity in the air hole size shows 
very little dependence of spectral response on the hole inhomogeneity. 

Polarization independence is reported to be a manifestation of rota- 
tional symmetry with no translational symmetry 16 . Orientation inde- 
pendence in the transmission enhancement is shown in Fig. 2(c) for 
different 0 and (p for TM polarized incident light. Measured launch 
angle dependent transmission enhancement spectra are plotted as a 
contour plot in Fig. 2(d) and 2(e) for TE and TM polarizations, 
respectively. Yellow contours corresponding to different SPP modes 
are identified from the calculations described above and are marked. 
ki mode is weak and the corresponding contours are not resolvable 
in the plots. Broad featureless red regions covering 500-700 nm 
wavelength range show the broadband transmission enhancement 
that is launch angle independent. 

In addition, the multiple periodicities and dense k-space displayed 
could lead to greater number of SPP modes in the same wavelength 
range compared to 1-d or 2-d periodic structures, resulting in broad 
band response. Enhanced transmission is known to happen due to 



coupling of field from top to bottom. Thus, in addition to the pattern, 
metal thickness is also critical for efficient coupling and transmission 
enhancement. To rule out the role of high hole density (and fill 
factor) in the observed transmission enhancement, we tested a struc- 
ture with 2-d square array of holes in metal with same metal thick- 
ness and with the same hole density (and thus fill factor) as the 
quasicrystal sample. Thus the square lattice has period of 150 nm 
with 85 ± 10 nm diameter of holes. Measured diffraction pattern 
(Fig. 3(a)) with 442 nm laser is much sparser and weaker compared 
to the quasicrystal pattern indicating that we may observe narrow 
and weak transmission enhancement in the case of square lattice. 
Accordingly, we observed much lower transmission enhancement of 
— 1.1 compared to unpatterned metal for this structure as shown in 
Fig. 3(b). Sparse k-space and dispersive modes are already reported 
in square lattice patterns with period around 0.7 um 38 . 

To check the role of inhomogeneity in the hole size, we compared 
two 1 um period quasicrystal structures with varying distribution in 
hole diameter. One structure has hole size distribution of 92.5 ± 
7.5 nm and the other has 115 ± 35 nm. Almost identical transmis- 
sion behavior but with reduced resolution could be seen in Figure 4 as 
expected when inhomogeneity increases. It may be noted that the dip 
at about 950 nm observed in the spectrum is related to Wood's 
anomaly as has been predicted in quasicrystals 40 . 

It is advantageous to have control over the design to realize PlQC 
structures with desired spectral response. This can have applications 
in light harvesting, filters, nonlinear optics and quantum plasmonics 
where manipulation of dipole emitters is needed. We show such 
feasibility in a dual grating PlQC realized by combining PlQCs with 
two different base periods of 0.6 and 1 urn in which the transmission 
band has been increased to 800 nm. To realize the bigrating, the 
1 urn base period QC pattern was superposed on 0.6 um base period 
QC pattern with a common center. Schematic of the individual pat- 
terns for the three QCs, 0.6 um, 1 um and the combined bigrating 
pattern, are shown in Fig. 5a, b, and c, respectively. Figure 5d shows 
the SEM image of the bigrating quasicrystal and Figure 5e shows the 
measured transmission enhancement curves for the three structures. 
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Up to 15 times enhancement in transmission over a broad wave- 
length range could be seen for the bigrating structure for both TE and 
TM polarizations. Thus, multi-grating PlQCs could be a way forward 
to realize structures with required transmission response. 

Conclusion 

We demonstrated a method to design and fabricate large area qua- 
sicrystal patterns of a given rotational symmetry by electron beam 
lithography. As an example, we demonstrated tt/5 rotation symmetry 
plasmonic quasicrystals with clear optical diffraction pattern. We 
show in PlQCs with two different periods that the inherent lattice 
constant (or period) in the present method allows one to design the 
spectral response. These structures showed near dispersionless 
enhanced transmission over a broad wavelength range for both TE 
and TM polarizations. Also, the present method can be extended to 
multi-grating patterns to achieve broadband spectral response. We 
demonstrated one such structure with very broadband transmission 
enhancement of over an order of magnitude in bigrating PlQC com- 
pared to unpatterned gold film. 

Methods 

To prepare the sample, a 100 nm thick gold film was deposited in Argon (Ar) plasma 
by radio frequency (RF) magnetron sputtering on optically flat 1 mm thick quartz 
substrate. A polymer, electron beam resist, PMMA 495 A4, was spin coated on top of 
the gold film at 2000 rpm to get a resist thickness of 200 nm. Electron beam litho- 
graphy (EBL) was then performed by using Raith e-line to pattern the resist with the 
designed quasicrystal pattern. A high density of holes meant that we had to use dose 
scaling to get good feature resolution. Subsequent development was done using 
MIBK:IPA (1:3) for 90 seconds followed by IPA treatment for 60 seconds. Reactive 
Ion Etching (RIE) was then used to Argon (Ar) sputter the exposed Gold regions to 
get air hole pattern in gold film. The Ar flow rate was 50 seem and chamber pressure 
was maintained at 1 Pa. The radio frequency (RF) power used was 138 W. Etching 
was done for 7 minutes to remove the gold. The resist was then stripped using 0 2 
plasma at 80 W RF power at a flow rate of 50 seem for 6 minutes. To etch metals we 
modified the Sentech SI 500 inductively coupled plasma (ICP) RIE by placing a Teflon 
shield to cover and protect the ICP head. We reduced the damage of the gold surface 
by having 1-2 minute etch stop for every 1 minute etching. 

For linear optical transmission, a 100 W tungsten halogen lamp was used as the 
white light source with multiple lenses and pinholes to collimate the beam to < 0.3° 
divergence. A fiber coupled spectrometer (either Ocean Optics S4000 with resolution 
0.26 nm or 0.5 m Acton SP2500 spectrometer with thermoelectrically cooled PIXIS 
400BR CCD with 0.05 nm resolution at 600 nm) was used to collect and characterize 
the signal in the 450-1200 nm wavelength range. A Glan-Thompson polarizer 
(extinction ratio 10 5 :1) was used to get linearly polarized light (TE or TM). To 
perform angle resolved measurements, a motorized rotation stage with 0.03° angle 
resolution was used. 

C-programs were written to create the design and extract the coordinates of the 
lattice points of the quasi crystal pattern as well as for Quicksort algorithm. 
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